Onset of mechanical stability in random packings of frictional spheres 
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Using sedimentation to obtain precisely controlled packings of noncohesive spheres, we find that 
the volume fraction 0rlp of the loosest mechanically stable packing is in an operational sense well 
defined by a limit process. This random loose packing volume fraction decreases with decreasing 
pressure p and increasing interparticle friction coefficient fi. Using X-ray tomography to correct for 
a container boundary effect that depends on particle size, we find for rough particles in the limit 
p ^ a new lower bound, (^rlp = 0.550 ± 0.001. 

PACS numbers; 83.80.Fg, 46.65. +g, 45.70.Cc, 47.57.ef 



Introduction. - If granular materials such as sand, 
sugar, or snow are excited strongly (e.g. by shaking or 
shearing), they exhibit fluid- like behavior. However, af- 
ter the excitation stops, dissipation quickly produces a 
static packing that is mechanically stable under its own 
weight. Experiments [iLi, i, H i, i, 0, il and simu- 
lations [13, [ll|i [3 [3 Hi I have shown that the volume 
fraction has a well defined lower limit, 0rlp7 called Ran- 
dom Loose Packing (RLP). 

The value of 0rlp depends on the particle-particle in- 
teractions. Packings of cohesive particles like fine pow- 
ders are stable under their own weight for values of (/irlp 
as low as 0.15 [1, However, many granular mate- 

rials do not exhibit cohesive forces. Simulations of fric- 
tionless elastic noncohesive spheres have found the onset 
of a finite bulk modulus at the jamming point, 0j k, 0.64 
[13, [lHj [i^] ■ Real spheres have friction and then it has 
been suggested that 0rlp depends on the density dif- 
ference between the particles and the surrounding fluid 
[U, 0, S, The lowest volume fraction reported thus 
far, (^RLP = 0.555, was observed for slowly sedimenting 
spheres in a liquid of nearly the same density Q . 

Here we demonstrate a limit process that yields well 
defined values of 0rlp that depend on pressure and coef- 
ficient of friction. The results are discussed in the context 
of a statistical mechanics approach based on the ensem- 
ble of all mechanically stable configurations (l5l |. 

Experiment. - Mechanically stable packings of glass 
spheres were prepared by allowing the particles to sedi- 
ment following flow pulses in a water fluidized bed. The 
fluidized bed was contained in a vertical polycarbonate 
tube with an inner diameter Z? of 12.8 mm and a length of 
230 mm. The tube's bottom end was closed by a distrib- 
utor consisting of a porous bronze disc (height, 8 mm; 
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FIG. 1: (Color online) The approach to Random Loose Pack- 
ing in a limit process is achieved using flow pulses in a liquid 
fluidized bed. Data for different particle diameters were fit to 
((T]) to obtain (fiRLP- Letters in parentheses refer to the parti- 
cle samples in Table HI Inset: Diameter dependence of (pni-p 
without container size correction (see text). Sample height 
was 97 mm at RLP; five flow pulses were averaged for each 
flow rate. 



nominal pore size, 25 fim). A programmable syringe 
pump (Harvard Apparatus) created pulses of constant 
flow rate Q. During a flow pulse of 2 minutes length the 
granular medium fluidized and expanded until it reached 
a steady state height. After each flow pulse, the particles 
sedimented onto the distributor and formed a mechani- 
cally stable packing whose volume fraction depended on 
Q, as shown in Fig. [TJ A higher value of Q resulted 
in a more expanded fluidized bed, longer sedimentation 
time, and lower (jj of the sediment. Packings created in 
this way are independent of the state of the sample prior 
to the last flow pulse [l6| . This property is important 
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TABLE L Properties of the difFerent samples of glass spheres. 
Particle diameters d and standard deviations g were measured 
with a Camsizer (Retsch Technology). Angles of repose under 
water were averaged over 10 measurements. 



sample 


d (j-im) 


a{%) 


supplier 


angle of repose 


A 


96 


15.6 


Cataphote 


24.8 ± 1.0 


B 


167 


16.1 


Cataphote 


26.1 ± 0.7 


C 


277 


7.6 


Cataphote 


25.3 ± 0.8 


D 


322 


9.3 


Cataphote 


25.5 ± 0.7 


E 


261 


5.0 


MoSci 


24.0 ± 0.8 


F 


257 


7.8 


Cataphote 


27.7 ± 1.3 


G 


255 


2.7 


MoSci 


26.6 ± 0.7 



for any statistical mechanics approach [l7l |. The volume 
fractions in Fig. [T] are averaged over the whole sample: 
(/) — m/ pAhscd, where p is the particle density, m is the 
total mass of all the spheres, A is the cross sectional area 
of the tube, and /igod is the height of the sedimented sam- 
ple determined from images. 

The properties of the different samples of particles are 
given in Table H] The density p of the particles in each 
sample was measured with an accuracy of 0.06 % using 
a Gay-Lussac specific gravity bottle and a Micromeritics 
AccuPyc 1330 gas pycnometer; the average p was 2.48 
g/cm"^. To characterize the frictional properties of the 
samples we measured the angle of repose under water: a 
beaker containing a layer of particles about 5 mm high 
was tilted until the particles started to move. To obtain 
an especially rough sample (F) we soaked part of sample 
C for 3 hours in hydrofluoric acid. Sample E consists of 
spheres that were smoothed by exposure to more than 
45,000 flow pulses in a fluidized bed p^ . 

RLP is defined by a limit process. - The main improve- 
ment over earlier studies using sedimenting particles [1,01 
is that our control of Q allows us to change the sedi- 
mentation time independent of the liquid density. This 
procedure reveals the convergence of (/) to (/irlp- The 
observation that the slowest relaxing preparation yields 
the loosest packin gs agree s with simulations of frictional 
discs and spheres Figure [T] shows that is 

well described by the fit function used in p^ . 

HQ) = 0RLP + (1) 

which we use to determine (/)rlp- 

Dependence on particle diameter. - The inset of Fig. [T] 
indicates that 0rlp decreases with particle diameter; 
however, this decrease is due to lower volume fraction 
near a container wall, an effect known since the earliest 
studies This effect is explained in Fig. [2ja): since 
particles cannot penetrate the container wall, voids are 
larger there and the volume fraction of the layer adjacent 
to the boundary is lower than (/)buik measured in the core 
of the sample. The difference between ^apparent averaged 
over the whole container and (j^huik increases with the ra- 
tio d/D and produces the trend displayed in the inset of 

Fig.m 



We examined the finite size effect using X-ray tomog- 
raphy [l^, Fig. [21 In each run positions of 1.5 x 10^ 
spheres were measured with a resolution of better than 
0.1% of a sphere diameter fl^. Figure [5]^c) shows the 
difference between the apparent volume fraction using 
all particles and the bulk volume fraction (measured us- 
ing the Voronoi volumes of all particles that are at 
least 4c? away from the container walls). For Q < 10 
ml/min, ^apparent valucs from the tomographic measure- 
ments agree with the results (dotted curve) for the 277 
fim particles in Fig.[Tl but for Q > 10 ml/min the tomo- 
graphic values are larger due to unavoidable vibrations 
during the recording of the tomogram. 

The inset of Fig.[2Ic) shows the difference between the 
bulk and apparent as a function of (/'apparent- A linear 
fit yields 

4>hn\k = '/'apparent + 0.122((/'apparont - 0.505). (2) 

For all further experiments we used only spheres with di- 
ameter 261 or 257 /im (sample E or F), and we corrected 
for the effect of finite container size using ^ . 

Influence of pressure. - The stress inside a column 
of grains differs from the hydrostatic case in two ways: 
(i) anisotropy - the horizontal stress CTxx in the column 
differs from the vertical stress a^z- (ii) wall friction - 
the part of the load carried by the frictional sidewalls 
increases with depth z below the surface. Consequently, 
tTxx and (Tzz saturate with z. In our analysis we use a pres- 
sure dependence on height given by the Janssen model 
[20| . which assumes a constant stress ratio K = a^^/azz 
everywhere in a sample. Experiments show that this 
model is a fair approximation in the absence of external 
loads [2l|, [l^l . The model gives a saturation of pressure 
(p = Czz) with depth, 

P(^)=ftat(l-e-^/'), (3) 

to a constant value Psat — ^pgD/AKp^jif, where /iw is the 
coefficient of particle-wall friction, g is the gravitational 
acceleration, Ap is the density difference between spheres 
and surrounding liquid, and I = D/AKp^^. Equation ^ 
indicates two ways of controlling the pressure distribution 
inside the column: 

I) Increasing the sample height, which increases the 
fraction of the sample at Psat • If (f'KLP increases with p, 
then the average 0rlp measured by our method should 
increase with sample height. This behavior is confirmed 
in Fig. O 

II) Decreasing the density difference Ap, which de- 
creases Psat but keeps the pressure profile unchanged. We 
increased the fluid density to as high as 2.39 g/cm'^, close 
to the 2.48 g/cm'^ particle density, by adding sodium 
polytungstate to the water. Results for different Ap 
(Fig. d]) again confirm that (/'rlp decreases with decreas- 
ing p. 

A limit 0|^Lp would be given by matching the fluid 
and particle densities, but in this limit there would be no 
sedimentation and no connected granular packings would 



3 




bulk particles 



boundary layer 



0.60 



0.59 



0.58 



0.57 - 



0.56 




10 




0.56 0.57 0.58 0.59 

'^'apparent 

bulk volume fraction 

apparent volume fraction 



20 30 40 

^\ow rate (ml/min) 



50 



FIG. 2: (a) Illustration of the necessity of a finite size correction of the volume fraction determined from the total sample volume. 
Particles near the boundary have a lower volume fraction than those in the bulk, (b) Cross section of a 3-dimensional x-ray 
tomogram of the fluidized bed; the particles have a diameter of 255 /xm (sample G). The inner diameter of the polycarbonate 
tube (black circle) is 12.8 mm. (c) Bulk volume fraction for particles that are at least 4ci away from boundary, and the apparent 



volume fraction 



."apparciit 



for all particles. The dotted line corresponds to the fit to the 277 fim particles in Fig. [T] The inset 



shows a linear fit ^ to the difi'erence between bulk and apparent volume fractions. Sample height was 39 mm at RLP. 
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FIG. 3: Dependence of on Q for different sample heights 
(261 fim spheres, sample E). Values of (p were corrected for 
finite size effect using (2). Solid lines are fits to (HJ. Inset: 
the resultant ^rlp values as a function of sample height. 



form. Therefore, we extrapolate to determine 4>^j^p- In 
the absence of theory we follow Q and use the pressure 
dependence close to the jamming point known for fric- 
tionless static soft spheres [l^, fTllll^. HJ] and frictionless 
thermal hard spheres [25j : 



^RLP 



(4) 



where we identify Ap ~ Psat- A fit of (/)rlp for smooth 
particles (sample E) in Fig.|4]yields (jiRhP — 

0.555±0.006. 

The value of a = 0.51 ± 0.25 is approximate because 
our derivation of dH) did not take into account the 
dependence of K [21| |. 
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FIG. 4: (^RLP as a function of the density difference between 
particles and fiuid, for smooth (261 fim) and rough (257 /im) 
particles. Solid lines are fits to (|4)|. Data have been corrected 
for the finite size effect using ((2)1. Sample height was 64 mm, 
and the points are averages obtained at the highest possible 
fiow rate, since for small density differences the fiow rate range 
was too small for a meaningful fit with ((l|. 



Influence of frictional properties. - Figure [4] shows 
that 0RLP for the rough spheres was lower than for the 
smooth spheres. For the rough spheres a fit to (|1]) yields 
(jjO^^p = 0.550 ± 0.001 and a = 0.89 ± 0.16. The decrease 
of ^RLP and </>j^Lp with increasing friction agrees with 
another experiment model [1^, and simulations [12, 

Discussion. - Our experimental results and numerical 
simulations [H, ll^l both show that RLP is well defined in 
an operational sense: in the limit of infinitesimally slow 
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preparation, the volume fraction of a sample converges 
to (f>RLP, independent of the details of preparation. Care 
should be taken in comparing theory for frictionless hard 
spheres with the experimental results, in part because of 
the different possible ways of defining mechanical stabil- 
ity [23, Hi. 

The observation of a well-defined (j)RLP can be con- 
sidered within the framework of a statistical mechanics 
of static granular material [l5[, where a configurational 
entropy S is defined as the logarithm of the number of 
mechanically stable configurations for a given (p, p, and 
friction coefRcient. Two different approaches can explain 
RLP using two different assumptions of how S depends 
on (f>. The first approach assumes that RLP is the small- 
est (j) where S becomes larger than zero. This is compat- 
ible with the existence of looser, highly ordered configu- 
rations [2?! , as their number seems not to grow exponen- 
tially with system size, so S* = 0. 

The second approach is supported by numerical results 
on the number of stable configurations of frictional discs 
[2^, where S has a maximum at RLP. This idea agrees 
with slow sedimentation leading to RLP: it is simply the 
most probable configuration. If the sedimentation speed 
is increased, the additional kinetic energy allows the sys- 



tem to explore the local energy landscape and find rarer 
but lower potential energy (denser) configurations. Fur- 
ther, the maximum of S and therefore RLP moves to 
higher values of for decreasing friction [23| . This agrees 
with our results and with simulations of frictionless disks 
that have a maximum of S at Random Close Packing 
[29'| . Our results indicate also that increasing p shifts the 
maximum of 5 in a similar way. 

Conclusions. - Mechanically stable packings of spheres 
prepared with increasing sedimentation time display a 
lower bound of their volume fraction, (/)rlp, which de- 
pends on the pressure and the coefficient of friction but 
not on the diameter of the spheres. In the limit of 
zero pressure we have found a new lowest value of 0rlp , 
0.550 ±0.001. 
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